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1. 


INTRODUCTION 


During  the  past  decade  a  number  of  costly  aircraft  accidents  have  occurred 
where  wind  shear  has  been  identified  as  the  probable  cause  .  The  identifica¬ 
tion  of  a  wind  shear  hazard  for  large  jet  aircraft  taking  off  or  on  final 
approach  has  prompted  the  Federal  Aviation  Administration  (FAA)  to  examine 
various  ground  based  in-situ  and  remote  sensing  devices  to  determine  if  a  wind 
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shear  warning  system  could  be  developed.  Surface  arrays  of  pressure,  temperature  , 
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and  wind  sensors  may  give  indications  of  dangerous  winds  near  the  airport ,  but 
remote  sensing  devices  offer  the  possibility  of  measuring  dangerous  winds  along 
the  actual  track  of  the  aircraft.  Early  research  work  to  solve  the  wind  shear 
problem  with  remote  sensing  led  to  the  installation  and  test  of  a  system  for 
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measuring  wind  profiles  at  Dulles  International  Airport  .  This  system  used  an 
acoustic  radar  for  measuring  vertical  profiles  of  horizontal  wind  to  500  m 
altitude  in  the  precipitation-free  atmosphere  and,  during  rain,  when  the  acoustic 
system  could  not  obtain  data  because  of  excessive  noise,  a  low-power  microwave 
radar  measured  wind  profiles.  However,  measurement  of  the  vertical  profile  of 
horizontal  wind  at  one  point  near  the  runways  can  only  detect  those  wind  shear 
hazards  that  arise  from  vertical  gradients  of  the  horizontal  wind.  Hazardous 
winds  from  thunderstorm  outflows,  for  example,  may  not  be  detected  by  such  a 
system.  Therefore,  research  work  has  continued  to  develop  remote  sensing  tech¬ 
niques  for  observing  winds  along  the  aircraft  track.  Both  laser  radar  systems5 
and  microwave  radar  systems  offer  this  potential. 

Microwave  radars  have  been  a  tool  for  operational  meteorology  for  over 
20  years.  During  the  1970's,  Doppler  radar  technology  has  developed  so  that 
the  next-generation  of  weather  radars  will  include  Doppler  capability6.  The 
Doppler  weather  radars  measure  the  reflectivity  of  the  hydrometeor  scatterers, 
just  as  conventional  or  non-Doppler  radars,  but  they  also  measure  the  average 
radi al  velocity  and  the  spread  of  the  radial  velocities  in  each  radar  resolution 
celJ .  If  the  Doppler  radar  is  sufficiently  sensitive,  it  can  measure  the  re- 
fle<  tivity  and  velocity  from  scattering  in  the  precipitation-free  atmosphere. 

In  t  his  case  the  scattering  arises  from  refractive  index  variations  which  occur 
threughout  the  lower  atmosphere  and  especially  in  the  planetary  boundary  layer. 
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The  radar  senses  that  component  of  refractive  turbulence  whose  scale  size  is 
one-half  the  radar  wavelength.  Typical  meteorological  radars  with  wavelengths 
of  5  to  20  cm  have  sufficient  sensitivity  to  detect  refractive  turbulence  in 
the  planetary  boundary  layer,  sometimes  to  a  range  of  more  than  50  km7.  Micro- 
wave  radar,  therefore,  offers  the  possibility  of  all-weather  wind  measurements 
to  altitudes  needed  to  insure  safe  operation  of  aircraft  in  the  terminal  area. 

The  actual  altitudes  where  wind  shear  poses  a  threat  to  aircraft  depends  on  the 
aircraft  response  and  therefore  on  the  type  of  aircraft.  Measurement  of  wind 
shear  in  the  lowest  500  m  has  been  the  goal  for  remote  sensing  systems. 

In  January,  1977  the  National  Oceanic  and  Atmospheric  Administration  (NOAA) 
and  the  FAA  began  a  joint  program  to  examine  pulse  Doppler  radar  as  a  potential 
wind  shear  detection  device.  [NOAA's  Wave  Propagation  Laboratory  (WPL)  had 
been  involved  with  the  Doppler  radar  work  at  Dulles  and  had  also  developed  a 
Frequency  Modulated- Continuous  Wave  (FM-CW)  Doppler  radar  for  wind  measurements 
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in  the  precipitation-free  atmosphere  .  The  FM-CW  radar  is  also  being  tested  as 
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a  potential  wind  shear  warning  device  for  the  United  States  Air  Force  ] .  The 
objective  of  the  joint  program  was  to  explore  whether  an  existing  FAA  air  traffic 
control  radar  could  be  configured  for  wind  measurements.  A  radar  wind  shear 
warning  system  would  be  much  more  attractive  if  existing  radars  could  be  used. 
Experimental  work  was  to  be  performed  at  the  National  Aviation  Facilities 
Experimental  Center  (NAFEC)  at  Atlantic  City,  NJ.  Limitations  of  Doppler  radar 
were  recognized  from  the  start  and  some  of  the  important  limitations  would  be 
explored  by  actual  testing  at  NAFEC  and  some  by  computer  simulations.  The  radar 
was  to  be  equipped  with  a  15  foot  parabolic  reflector  antenna  and  a  simple  data 
processing  system,  sufficient  to  demonstrate  the  ability  of  radar  to  measure 
Doppler  spectra  in  clear  air,  would  be  used.  Major  questions  to  be  answered 
with  the  radar  were:  (a)  how  would  ground  clutter  restrict  the  minimum  range 
or  altitude  at  which  clear  air  radar  signals  could  be  measured  and  (b)  was  the 
radar  sensitive  enough  to  measure  clear  air  returns  under  all  conditions.  Major 
questions  to  be  answered  by  computer  simulation  were:  (a)  does  a  15  foot 
diameter  antenna  give  sufficient  resolution  to  measure  wind  profiles  and  (b)  how 
should  the  radar,  which  measures  only  one  component  of  the  wind,  be  used  to 
detect  and  measure  dangerous  winds. 
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2. 


NAFEC  RADAR  SYSTEM  AND  TESTS 


One  channel  of  an  ASR-8  radar  has  been  configured  for  wind  shear  tests. 

Figure  2.1  shows  a  block  diagram  of  the  system.  The  15  foot  parabolic  re¬ 
flector  was  procured  from  the  Microwave  Specialty  Corp.  in  San  Diego,  CA  and 
was  installed  in  the  spring  of  1978.  A  pedestal  was  not  available  until  May, 

1979  so  many  of  the  tests  were  conducted  with  a  fixed  antenna  position.  The  trans¬ 
mitter,  preamplifier  and  mixer  of  the  ASR-8  were  used  without  modification.  Al¬ 
though  the  ASR-8  has  a  linear  IF  and  in-phase  and  quadrature  phase  detectors,  they 
were  not  used  because  a  better  receiver  chain  was  available  from  a  previous  NAFEC 
project.  WPL  provided  the  range  gate  generator  and  buffer  memory  on  a  NOVA  com¬ 
puter  board.  The  computer  and  its  peripherals  were  available  from  previous  pro¬ 
jects.  The  IMLAC  graphics  display  proved  unreliable  and  was  replaced  with  a 
Tektronix  graphics  terminal.  Software  to  observe  the  Doppler  spectra  at  8  range 
locations  was  developed.  The  Doppler  spectra  are  calculated  in  software  and  are 
measured  for  only  one  range  interval  at  a  time.  This  was  recognized  as  a  limita¬ 
tion  but  it  was  deemed  suitable  for  demonstration.  [Equipment  for  calculating  the 
spectra  for  16  or  more  range  intervals  simultaneously  is  commercially  available. 
Although  the  mean  velocity  could  be  calculated  for  all  range  intervals  using  the 
so-called  pulse-pair  technique,  the  complete  Doppler  spectra  is  needed  for  these 
tests  where  ground  clutter  is  expected  to  be  severe.]  The  software  enables  the 
operator  to  select  the  range  delay  to  the  first  range  interval,  the  spacing  between 
intervals,  the  number  of  range  locations  (<  8),  the  number  of  spectra  that  would  be 
averaged,  and  the  scaling  and  baseline  of  the  display.  The  number  of  points  '.n  the 
Doppler  spectrum  was  fixed  at  128.  The  unambiguous  Doppler  velocity  is  ±  27.55  m/s 
or  about  0.43  m/s  for  each  spectral  point.  The  first  tests  of  the  system  were  made 
in  June,  1978  when  the  radar-measured  radial  velocity  profiles  were  compared  with 
measurements  made  with  an  aircraft1®.  A  second  series  of  tests  comparing  radar  and 
aircraft  measurements  has  recently  been  made11.  These  tests  demonstrate  that  the 
radial  component  of  wind  can  be  measured  by  Doppler  radar  in  various  types  of 
weather.  Not  enough  tests  have  been  made  to  determine  if  measurements  can  be  made 
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in  all  weather  regimes.  However,  measurements  with  a  less  sensitive  radar  in¬ 
dicate  that  the  NAFEC  system  should  be  able  to  measure  radial  wind  profiles  in  the 
boundary  layer  practically  all  the  time. 


Figure  2.1.  Block  diagram  of  NAFEC  wind  shear  detection  system. 
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RADAR  MEASUREMENT  OF  WIND  PROFILES 


If  a  radar  is  positioned  at  the  end  of  a  runway  and  is  pointing  along  a 
3  deg  glide  slope,  it  should  be  able  to  measure  the  headwinds  encountered  by 
the  aircraft.  (Headwind  will  be  used  to  denote  a  headwind  or  tailwind;  negative 
headwinds  are  tailwinds.)  However,  if  there  is  vertical  shear  of  this  wind  com¬ 
ponent,  the  radar  will  measure  the  actual  profile  only  if  its  resolution  is  suffi¬ 
cient.  For  example,  an  antenna  with  a  one-way  beamwidth  of  1.5  deg,  (a  value  similar 
to  that  of  the  NAFEC  test  radar)  would  have  a  2-way  linear  beamwidth  of  (0.018) R 
where  R  is  the  range.  At  a  3  deg  elevation  angle,  the  2-way  linear  beamwidth 
is  0.342  h  where  h  is  the  height.  The  actual  wind  profile  will  be  filtered 
or  smoothed  by  the  radar  resolution  volume. 

In  phase  1  of  this  contract  computer  simulations  were  made  to  examine  if  a 

radar  with  a  1.5  degree  beamwidth  could  resolve  wind  profiles  used  to  represent 
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hazards  to  low-level  aircraft  .  Figure  3.1  shows  typical  results.  The  actual 
wind,  shown  by  solid  lines,  varies  only  with  height.  The  expected  value  of  the 
radar-measured  wind  profile  is  shown  with  diamond  points.  Note  that  the  (unrea¬ 
listic)  gradients  of  the  input  profile  are  filtered  by  the  vertical  dimensions  of 
the  radar  beamwidth.  When  the  elevation  angle  is  increased  from  3  deg  (a)  to  6  deg 
(b) ,  the  vertical  beamwidth  at  each  measurement  height  is  reduced  by  2  so  the  radar 
resolution  improves.  The  resolution  continues  to  improve  as  elevation  angle  is 
increased  until  the  vertical  extent  of  the  resolution  along  the  beam  becomes  im¬ 
portant.  A  pulse  length  of  1  microsecond  (range  resolution  of  150  m)  was  used  in 
the  simulations.  The  radial  resolution  becomes  more  important  them  the  cross-beam 
resolution  for  heights  of  500  m  or  less  when  the  elevation  angle  is  greater  than  14 
deg.  Therefore  the  ability  of  the  radar  to  resolve  low-level  wind  profiles  will  be 
best  for  elevation  angles  somewhat  less  than  14  deg. 

The  results  from  the  simulations  show  that  the  radar  can  adequately  resolve 
low-level  wind  profiles  of  the  type  that  can  be  hazardous  to  aircraft.  The 
Doppler  radar  can  therefore  serve  the  following  functions: 

a.  measure  vertical  profiles  of  horizontal  wind  to  500  m  altitude  or  more 
whenever  the  wind  field  varies  only  with  altitude 
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b.  measure  vertical  profiles  of  horizontal  wind  to  500  m  altitude  or  more 
whenever  the  wind  field  varies  with  altitude  and  varies  linearly  with 
horizontal  distance 

c.  measure  the  headwinds  along  approach  and  depart  lire  paths  whenever  the  radar 
is  located  such  that  the  aircraft  track  is  approximately  toward  or  away 
from  the  radar.  The  headwind  component  is  more  critical  than  the  crosswind 
component  for  landing  aircraft  because  changes  in  headwind  cause  air  speed 
changes  which  cause  approaches  to  be  long  or  short 

d.  if  the  radar  is  not  suitably  located  for  measuring  headwinds  along  the 
aircraft  track,  the  radar  can  measure  the  radial  velocity  profile  along  the 
approach  path  and  check  to  see  if  it  is  consistent  with  the  vertical  profile 
of  horizontal  wind. 

A  single  Doppler  radar  can  measure  vertical  profiles  of  horizontal  winds  for 
uniform  or  linear  conditions;  however,  the  profile  measured  at  the  radar  site 
is  not  applicable  to  the  aircraft  track  located  some  distance  away  if  there 
are  horizontal  gradients.  If  the  horizontal  wind  field  varies  only  with 
height,  then  the  radar-measured  wind  profile  can  be  directly  applied  to  land¬ 
ing  aircraft.  Therefore,  the  radar  can  be  used  to  measure  wind  profiles  as  in 
(a)  and  (b)  and  to  check  the  applicability  of  the  measured  profile  as  in  (c) 
or  (d) . 

The  most  widely  used  method  of  measuring  wind  profiles  with  Doppler  radar 
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is  the  velocity-azimuth  display  (VAD)  method  .  The  radial  velocity  is  mea¬ 
sured  at  each  height  as  a  function  of  azimuth  (continuously  or  discretely)  with 
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a  fixed  antenna  elevation  angle.  Fourier  analysis  of  the  radial  velocity 

measurements  gives  not  only  the  mean  wind,  but  also  the  shearing  and  stretching 

3u  3v  du  3 

deformations,  +  -^)  and  -  -gy-) ,  and  the  vertical  motion  and/or  divergence 

{■jj^  +  .  Here,  u,  v,  w  are  the  wind  components,  which  will  be  averages  over 

a  circle  of  radius  z/tan8  where  z  is  the  measurement  height  and  9  is  the  antenna 

E  E 

elevation  angle. 
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A  simplified  version  of  the  VAD  scan,  used  in  the  Dulles  wind  shear  system  , 
is  useful  for  measuring  mean  wind  profiles  routinely.  This  method  involves  mea¬ 
suring  the  radial  velocity  profile  at  4  azimuths  separated  by  90  deg.  This  method 
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allows  the  use  of  fixed  antennas  which  may  be  less  costly  than  a  single  antenna 
and  pedestal.  Let  the  tu,  v]  wind  components  be  defined  by  the  orthogonal  pointing 
directions.  Then  the  4  radial  velocity  measurements  are: 


(h)  =  u^h)  cos6e 

V2  (h)  =  v2  (h)  cos0E 

V3(h)  =  -u3(h)  cos0E 

V.  (h)  =  -v  (h)  cos0 
4  4  E 


+  (w  +  Vt)1(h)  sin©E 

+  (w  +  Vfc)2(h)  sin6E 

+  (w  +  vfc) 3 (h)  sin©E 

+  (w  +  vt)4(b)  sin0E. 


is  the  fall  speed  of  the  scatterers  when  there  is  no  vertical  wind.  V  is 
zero  for  refractive  turbulence  scattering  but  can  be  as  large  as  9  m/s  in  pre¬ 
cipitation15  and  even  larger  if  there  are  hail  particles. 


If  the  vertical  motion  is  assumed  to  be  the  same  at  the  4  measurement  points 
(or  if  the  elevation  angle  is  low) ,  then  the  u  and  v  profiles  are  found  from: 


u(h) 

v  (h) 


u^  (h)  +  u3(h)  V1(h)  -  V3(h) 

2  2  cos0 

£ 

v2 (h)  +  V4 (h)  V2 (h)  -  V4 (h) 

2  "*  2  cos0_ 

£ 


Thus  under  this  assumption  the  wind  profile  is  measured  correctly  for  a  uniform 
horizontal  wind  or  a  horizontal  wind  that  varies  linearly  with  distance.  The 
error  in  u  or  v  caused  by  changes  in  vertical  motion  between  opposite  measurement 
points  is  [A(w  +  Vt)][tan0E]  or  (h) [gradient (w  +  V^)].  The  vertical  motion  cannot 
be  found  unless  the  horizontal  wind  field  is  uniform.  Just  as  in  the  complete  VAD 
analysis  technique,  the  vertical  motions  cannot  be  separated  from  gradients  of 
horizontal  wind.  This  simplified  method  of  wind  profile  measurement  is  preferred 
to  the  full  VAD  technique  for  precipitation-free  conditions  because  of  the  long 
dwell  times  needed  at  each  azimuth.  The  largest  error  of  horizontal  wind  caused  by 
the  fall  speed  of  raindrops  occurs  if  the  scattering  in  one  direction  is  from 
refractive  turbulence  or  very  small  raindrops  with  near-zero  fall  speed  and  the 
scattering  in  the  opposite  direction  is  from  large  raindrops.  The  maximum  error 
would  b»  9  tan0E  or  less  than  2  m/s  for  elevation  angles  as  high  as  12  deg.  Note 


the  sum  of  the  radial  velocity  measurements  made  in  opposite  directions  will 
be  a  measure  of  vertical  particle  velocity  and  horizontal  wind  gradient.  That 
is: 


V1(h)  +  V3<h)  =  [ux(h)  -  u3(h)]  cos0E  +  [(w  +  Vfc)1  +  (w  +  Vfc)3]  sin6E  . 

Since  the  contributions  to  this  sum  from  the  fall  speed  of  raindrops  and  normal 
vertical  air  motions  are  bounded,  this  sum  should  not  deviate  from  zero  by  more 
than  a  few  m/s.  If  it  does  it  is  an  indication  of  extreme  vertical  motions 
such  as  from  hailstones  or  strong  vertical  drafts  or  an  indication  of  horizontal 
wind  gradients.  The  sum  of  the  radial  measurements  is  one  check  of  whether  the 
measured  profile  is  applicable  to  points  other  them  a  line  above  the  radar. 
Additional  tests  of  whether  the  measured  profile  is  representative  of  the  wind 
sente  distance  away  from  the  radar  can  be  made  by  measuring  the  wind  profile  at 
different  elevation  angles  and  different  azimuth  angles. 

Another  check  of  whether  the  radar-measured  wind  profile  is  valid  along  the 
aircraft  track  can  be  made  by  measuring  the  radial  velocity  along  the  track, 
and  comparing  the  measurements  with  the  expected  value.  That  is,  the  measured 
velocity  along  the  aircraft  track  is  u  cosfj  cos0„  +  v  sin3  cos0  +  (w  +  V  )  sin0  . 

E  E  t  E 

The  measurement  must  be  consistent  with  the  radial  velocity  predicted  with  this 

equation  with  u,  v,  and  (w  +  V  )  sin0  found  from  the  4-point  profile  measure- 

u  E 

ment. 


The  conclusions  from  phase  1  simulations  were  that  radar  wind  profiles 
could  be  measured  by  Doppler  radar  and  that  a  radar  located  near  the  inter¬ 
section  of  major  runways  could  measure  the  profile  of  headwind  along  approach 
and  departure  paths.  In  general,  however,  a  single  radar  is  not  able  to  measure 
winds  along  the  approach  or  departure  path  of  all  runways  and  wind  profiles 
measured  at  the  radar  cannot  be  applied  to  the  aircraft  track  in  the  case  of 
hazardous  winds  from  localized,  transient  phenomena. 
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4. 


RADAR  DETECTION  OF  WIND  SHEAR 


A  Doppler  radar  can  measure  only  the  radial  component  of  the  wind  so  it 
cannot,  in  general,  measure  the  horizontal  winds  along  the  aircraft  track. 

If  one  radar  is  to  provide  wind  shear  warning  data  for  the  entire  airport  it 
must  be  able  to  detect  hazardous  winds  regardless  of  its  location  relative 
to  any  one  approach  path.  Transient  events,  such  as  gust  fronts  or  outflow 
from  thunderstorms  produce  locally  complex  wind  fields  that  cannot  be  directly 
measured  with  a  single  Doppler  device. 

Complex  three-dimensional  wind  fields  can  be  measured  with  suitably 

placed  multiple  Doppler  radars.  However,  the  data  processing  used  to  achieve 

these  wind  fields  is  performed  off-line  and  multiple  Doppler  systems  are 

generally  viewed  as  research  devices.  One  of  the  research  outputs  of  multiple 

Doppler  experiments  (one  that  has  been  largely  ignored)  is  to  determine  what 

information  about  the  complex  dynamics  of  atmospheric  motion  could  be  derived 

from  a  single  Doppler  radar.  When  (single)  Doppler  radars  become  operational 

tools,  us  expected  in  the  next  generation  of  weather  radars,  there  should  be 

increased  emphasis  on  deriving  maximum  utility  from  single  radar.  While  a 

complete  dual-Doppler  radar  system  which  derives  horizontal  wind  fields  in 

real  time  could  be  constructed,  (using  two  radars  located  10-20  km  from  the 

airport  such  that  2-dimensional  data  for  all  approach  and  departure  paths 

would  be  obtained)  such  a  system  would  be  complex  and  would  have  relatively 

poor  altitude  resolution  unless  the  antenna  apertures  were  very  large.  A 

possible  solution  for  the  wind  shear  problem  would  be  to  locate  radars  near 

the  intersection  of  major  runways  so  that  the  headwind  component  could  always 

be  observed.  The  radars  could  use  a  common  data  system  so  only  the  transmit/ 

receive  units  and  the  antennas  would  need  to  be  at  the  runways.  If  these 

runway  cadars  were  required  to  operate  only  during  precipitation ,  then  low 

power  and  inexpensive  radars  (X-band  with  aim  diameter  antenna)  for  each 

runway  '-ould  be  a  practical  solution.  These  radars  would  complement  a  single 

scanning  radar  that  could  measure  wind  profiles  in  all-weather.  If  the  runway 

radars  '/ere  required  to  measure  winds  in  all  weather  situations,  the  FM-CW  radar 

g 

may  provide  a  practical  solution  .  In  this  case,  the  runway  radars  would  complement 
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a  single  scanning  Doppler  radar  that  would  not  have  to  measure  wind  profiles 
in  all-weather.  Rather,  it  could  be  a  NEXRAD^  radar,  located  some  distance 
from  the  airport,  and  used  primarily  for  storm  surveillance . 

Efforts  have  been  made  to  use  data  from  a  single  radar  to  obtain  two 

dimensional  wind  fields.  One  method,  tested  on  a  convective  storm,  measured 

the  motion  of  the  reflectivity  patterns  by  correlating  reflectivity  fields 

measured  about  2  minutes  apart^®.  The  motion  of  fine  scale  patterns  in  the 

radial  velocity  field  can  also  be  measured  by  correlation  of  radial  velocity 

17 

fields  measured  at  different  times  .  Experiments  in  the  optically  clear 
air  show  that  these  eddy  fields  remain  correlated  so  that  their  motion  can  be 
measured.  The  spatial  resolution  for  winds  derived  from  correlation  methods 
will  be  much  poorer  them  for  dual  Doppler  measurements  because  relatively 
large  sections  of  data  must  be  used  to  obtain  the  motion  of  a  feature  or  per¬ 
turbation  of  the  radial  velocity  field  or  reflectivity  field.  Other  methods 

of  obtaining  wind  fields  from  single  Doppler  radar  involve  combining  data 

18  19 

obtained  from  a  sector  scan  or  volume  scan  .  Since  these  methods  also 
use  data  from  many  radar  resolution  volumes  to  derive  wind  field  data,  the 
resolution  will  be  relatively  poor  and  not  suitable  for  low-level  wind  shear 
detection . 

If  a  single  radar  measures  the  radial  velocity  at  2  points  near  the  glide 
slope,  then  a  2-dimensional  wind  vector  can  be  calculated  if  vertical  motion 
is  neglected.  A  radar  located  relative  to  the  runway  approach  path  as  in  Pig.  4.1 
cannot  directly  measure  the  headwind  component  encountered  by  the  aircraft.  How¬ 
ever,  by  measuring  the  radial  velocity  at  A  and  B,  the  longitudinal  velocity  (V  ) 
and  transverse  velocity  (VT)  can  be  derived  at  P  so  that  both  headwind  and  cross- 
Wind  can  be  found  on  the  aircraft  approach  path.  At  first  glance  this  appears  to 
be  a  method  which  does  not  rely  on  assuming  a  uniform  wind  over  a  large  distance, 
but  with  further  analysis  it  is  seen  that  the  method  fails  if  there  are  wind  gra¬ 
dients.  First,  assume  that  the  winds  are  locally  constant;  that  is,  the  wind 
is  the  same  at  A,  P,  and  B.  Let  the  antenna  elevation  angle  be  zero  to  simplify 
the  equations.  The  transverse  velocity  at  P  is  given  by 


Aircraft  Landing  Path 


Figure  4.1.  Single  Doppler  radar  measurement  of 
two-dimensional  winds. 
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where  V  and  V  are  the  radar -measured  radial  velocities  at  A  and  B.  The  longi- 
A  B 

tudinal  velocity  at  P  is  given  by 


VA  +  VB 
2  cos(-^p-) 


or  it  can  be  measured  directly  at  P  by  the  radar .  The  uncertainty  in  measurement 
of  VT  is  given  by 


Av_ 


AV„ 


sin(~) 


where  AV  is  the  uncertainty  in  the  measurement  of  radial  velocity  which  is 

K 

assumed  to  be  the  same  at  A  and  B.  Thus,  if  the  radial  velocities  are  measured 
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with  an  uncertainty  of  0.5  m/s,  the  transverse  velocity  is  measured  with  an  un¬ 
certainty  of  2  m/s  for  AAZ  =  20  deg.  If  AAZ  is  smaller,  then  the  assumption  of 
locally  constant  wind  is  made  for  a  smaller  area,  but  the  error  in  transverse 
velocity  increases  for  a  given  measurement  accuracy. 


However,  suppose  there  are  wind  gradients  and  the  two-point  method  is  used 
because  it  is  believed  that  if  the  measurements  are  made  at  closely  separated 
points  then  the  locally  constant  assumption  would  still  be  valid.  The  transverse 
velocity  is  given  by 


VT  = 


VA- 


B 


•  /AAZ, 
sin  (— 2 — ) 


2  k 


LT 


,AAZ, 
cos (-£— ) 


where  R  is  the  distance  from  the  radar  to  P  and  k  is  the  average  change  in  V 

LT  L 

per  unit  transverse  distance.  In  the  limit  of  small  angles. 


V 


T 


VA“ 


V„ 


Aaz 


+  2  k 


LT 
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so  that  no  matter  how  small  the  measurement  separation  the  gradient  term  introduces 
errors.  Decreasing  the  separation  angle  not  only  increases  the  measurement  un¬ 
certainty  (first  term)  but  also  increases  the  errors  caused  by  gradients.  R  is 
of  the  order  10  to  10  m  and  k^T  is  10  to  10  s  in  dangerous  shear  so  the 
errors  are  too  large  for  useful  measurements  whenever  there  are  dangerous  wind 
shears.  Thus,  a  single  Doppler  radar  cannot  adequately  measure  the  winds  en¬ 
countered  by  the  aircraft  in  those  cases  where  complex  3-dimensional  wind  fields 
are  expected. 


Even  if  a  single-Doppler  radar  cannot  measure  the  wind  profile  along 
the  aircraft  track,  it  may  be  able  to  detect  dangerous  winds  from  radial 
velocity  measurements.  The  profile  of  radial  wind  along  the  aircraft  track 
and  the  fie  Id  of  radial  velocity  measurements  in  the  space  near  the  track  may 

be  sufficient  to  detect  dangerous  winds.  Radar  measurements  provide  V  and  its 

>V  3v  9V  r 

it  r  r 

gradients  (  jp-,  -j-j-p,  -jj^)  •  If  the  radial  velocity  fields  measured  by  a  scanning 

E 

Doppler  radar  are  displayed  in  color  on  a  PPI  indicator,  characteristic  features  of 

j 
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certain  types  of  velocity  fields  are  readily  identified  .  Contours  of  constant 
radial  velocity  are  radial  lines  (Fig.  4.2)  if  the  wind  field  is  uniform. 
Non-uniform  wind  fields  cause  a  distortion  of  these  contours.  The  color  dis¬ 
plays  showing  mesocyclone  circulation  in  torhadic  storms6  are  the  best  example 
of  the  detection  of  hazardous  winds  by  single  Doppler  radar.  Another  display 

technique  that  has  proven  useful  for  single  Doppler  radar  data  is  the  multi- 

21 

moment  vector  display  .  Features  such  as  outflow  from  symmetrical  down- 

drafts  or  "downbursts"  should  be  easily  detected  in  the  pattern  of  the  radial 
12 

velocity  field  .  If  a  single  Doppler  radar  is  to  provide  a  wind  shear  warning 
system  for  an  entire  airport  runway  complex  (where  it  cannot  point  along  all 


Figure  4.2.  Contours  of  constant  radial  velocity  for  a  uniform 
wind  field  as  they  would  appear  on  a  PPI. 


runway  and  departure  paths) ,  then  it  must  be  able  to  detect  dangerous  wind 
gradients  in  the  radial  velocity  field.  A  computer  simulation  study  has  been 
made  to  determine  if  the  radar-measured  radial  velocity  could  detect  wind 
shear  in  wind  field  models  used  to  simulate  low-level  wind  shear  hazards. 


5 .  SIMULATION  PROGRAM 


12 

The  computer  program  used  in  phase  1  simulations  was  expanded  to  calculate 
and  display  the  actual  crosswind  and  vertical  wind  (in  addition  to  the  head¬ 
wind)  encountered  by  the  aircraft  on  the  approach  path.  The  expected  value 
of  the  radar-measured  radial  velocity  is  also  calculated  and  displayed.  In 
phase  1  simulations  the  radar-measured  radial  velocity  was  calculated  by 
integrating  the  actual  point  radial  velocity  function  over  the  entire  radar 
resolution  volume.  The  purpose  of  this  integration  was  to  test  if  the  radar 
could  resolve  the  existing  wind  profile.  Although  the  integration  remains 
part  of  the  simulation  program,  it  was  generally  not  used  for  these  tests. 

Rather,  the  radar-measured  radial  velocity  was  usually  taken  as  the  point 
value  at  beam  center.  Phase  1  simulations  showed  that  realistic  low-level 
profiles  could  be  resolved;  these  simulations  cure  intended  to  determine  if 
wind  shear  can  be  detected  in  the  radial  velocity  field. 

Figure  5.1  shows  the  geometry  for  the  simulation  program .  (The  setup 

is  similar  to  that  used  for  phase  1  simulations.  In  the  present  case  the 
location  and  orientations  of  the  runway  are  also  variable.)  Radar  range, 
azimuth,  and  elevation  (R,  AZ,  EL)  are  related  to  the  (x,  y,  z)  coordinates 
as  follows: 


x  =  R  sinAZ  cosEL  +  X„ 

R 

y  =  R  COSAZ  cosEL  +  Yr 
z  =  R  sinEL  +  ZR 

where  (X  ,  Y  ,  z  )  are  the  coordinates  of  the  radar  location.  The  aircraft 
K  R  R 

position  (X  ,  Y  ,  Z  )  is  given  by: 

G  G  G 
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mmm 


XG  =  RC  C°S<tG  Sin0G  +  XT 
Yg  =  Rg  cos<J>G  cos6g  +  Yt 

h  "  ZG  =  RG  sln*G  +  ZT 

where  (XT,  Y^,  Z^)  are  the  coordinates  of  the  touchdown  point,  RQ  is  the  dis¬ 
tance  from  the  aircraft  to  touchdown,  <J>_  is  the  glide  slope  angle,  0  is 

G  G 

the  runway  orientation,  and  h  is  the  height  of  the  aircraft.  The  runway 

orientation,  0_,  is  {10  x  runway  number  -  180)  deg.  The  program  allows  for 

arbitrary  location  of  both  the  radar  and  the  runway  relative  to  a  (u,  v,  w) 

wind  field  model  that  is  referenced  to  the  origin.  In  all  the  simulations  the 

glide  slope  angle  was  3  deg,  and  (X  ,  Y  ,  Z  )  and  Z  were  taKen  as  zero. 

R  R  R  T 

The  headwind,  crosswind,  and  vertical  wind  encountered  by  the  radar 

(V  ,  V  ,  W)  are: 

G  C 

VG  =  uQ  sin0G  cos <J>G  +  vG  COS0G  cos^  +  wQ  sin<J>G 

V  =  v  sin0  -  u  cos0 
c  G  G  G  G 

W  -wQ, 

where,  in  this  notation  u  ,  v  ,  w_  are  the  u,  v,  w  wind  components  on  the 

G  G  G 

glide  slbpe;  i.e.  ,  u  =  u(X  ,  Y  ,  z  )  etc. 

G  G  G  G 

The  position  of  the  aircraft  in  radar  coordinates  is 


»A  ■  ‘“c  '  V2  ♦  «G  -  IR>2  +  UG  -  ZR,2)V2 


A  =  tan"1  [ 


XG  ~  XR. 
YG  '  VR 


sin  ^  [ 


Z  -  Z 
G  R 


■]  . 
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The  radial  wind  component,  that  component  observed  by  the  radar,  is 


V  =  u  sinAZ  cosEL  +  v  cosAZ  cosEL  +  w  sinEL  , 

R 

where  (u,  v,  w)  are  the  wind  components  at  (R,  AZ,  EL) .  In  particular,  along 
the  glide  slope,  the  radial  velocity  is 

u_  sinA  cosE  +  v_  cosA  sinE  +  w_  sinE. 

G  G  G 

The  input  wind  models  are  in  Cartesian  coordinates  and  the  wind  components  and 
radial  wind  at  any  point  are  calculated  using  the  above  equations. 

One  of  the  computer  outputs  is  a  plot  of  the  actual  wind  (V  ,  V  ,  W)  and 

the  radial  wind  component  along  the  glide  slope.  The  simulation  program  also 
calculates  the  radial  velocity  field  as  would  be  observed  by  a  scanning  Doppler 
radar,  in  all  cases,  the  radar-measured  velocity  was  taken  as  the  radial 
velocity  at  beam  center  to  avoid  excessive  calculation  time.  Two  types  of 
scans  were  programmed:  a  volume  scam  and  a  constant  elevation  scan.  Ibe 
volume  scan  was  centered  at  the  point  where  the  aircraft  position  is  250  m 
above  touchdown.  The  scan  covers  an  azimuth  and  elevation  sector  to  scam  the 
entire  glide  path  for  h  =  0  to  500  m  and  covers  12  km  in  range.  The  azimuth 
sector  size  depends  on  the  radar  location  and  is  typically  60  deg.  The  com¬ 
puter  selects  the  required  azimuth  sector  (from  40  to  180  deg  in  20  deg  incre¬ 
ments)  .  The  elevation  scans  extend  to  the  elevation  angle  required  to  inter¬ 
sect  the  500  m  height  of  the  glide  slope.  The  scan  increments  are  0.1  deg 

in  elevation,  300  m  in  range,  and  the  azimuth  sector  is  covered  with  40 
increments.  The  output  can  be  displayed  as:  (a)  contours  of  constant  radial 
velocity  as  a  function  of  azimuth  and  the  height  that  the  seem  intersects 
the  glide  slope,  or  (b)  contours  of  constant  radial  velocity  as  a  function  of 
range  and  the  intersection  height.  Small  arrows  on  these  plots  indicate  where 
the  glide  slope  intersects  the  data  surface. 

The  constant  elevation  scan  is  also  centered  at  the  point  where  the  air¬ 
craft  position  is  250  m  above  touchdown.  Here,  the  scan  intersects  the  glide 
slope  at  only  1  point.  The  output  is  displayed  as  contours  of  constant  radial 
velocity  as  a  function  of  range  and  azimuth. 

18 


Realistic  runway  configurations  and  radar  locations  such  as 
Fig.  5.2  (representing  O' Hare  International  Airport)  were  used. 
4R  and  9L  were  used  as  representative  of  approach  paths  relative 
location.  The  simulation  program  is  described  in  Appendix  I. 


shown  in 
Runways  14R, 
to  the  ASR-8 


6 .  WIND  MODELS 

Wind  models  available  in  the  computer  simulation  program  include  both 
analytic  and  numerical  models  of  1,  2,  and  3  dimensional  wind  fields.  Models  1- 
6  are  the  same  models  used  in  phase  1  simulations. 

Models  1-3  are  1-dimensional  models  or  vertical  profiles  of  the  horizontal 

wind.  The  vertical  wind  is  zero.  These  models  were  developed  by  Fichtl  and 
22 

Camp“  .  Model  1  is  a  logarithmic  profile  representing  the  wind  in  the  neutral 
boundary  layer.  Model  2  represents  frontal  winds  and  model  3  represents  the 
nighttime  stable  boundary  layer. 

Model  4  (also  from  Fichtl  and  Camp)  represents  the  wind  profile  of  cold 

air  outflow  from  a  thunderstorm  along  the  glide  slope.  It  includes  vertical 

winds  but  it  is  not  a  thunderstorm  model;  it  represents  the  winds  along  a 

23 

line  through  a  storm.  Model  5  is  a  vertical  profile  of  horizontal  wind 
that  is  constant  in  direction  and  with  no  vertical  wind.  It  was  found  to  be 
difficult  for  pilots  to  handle  in  flight  trainers.  Model  6  is  a  3-dimensional 
analytic  model  of  an  axially  symmetric  downburst  from  a  thunderstorm.  The 
center  of  the  downburst  (X^ ,  Yc)  can  be  positioned  at  any  point  relative  to 
the  origin. 

Models  7-12  are  variations  of  model  6.  Model  12  has  the  same  functional 
form  as  model  6  but  with  different  scale  height  and  wind  speed.  Model  7  is 
the  sum  of  models  1  and  12.  Models  8-11  are  model  12  with  constants  added 
to  the  u  and  v  components,  i.e.,  am  axially  symmetric  downburst  with  transla¬ 
tion.  The  downburst  is  translated  by  u  =  10  m/s  in  model  8,  u  =  20  m/s  in 
model  9,  v  «*  10  m/s  in  model  10  and  v  =  20  m/s  in  model  11.  (The  translation 
is  constant  with  height  above  10  m  altitude  and  below  10  m  it  decreases  linearly 
to  zero  at  the  surface. 
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Figure  5.2.  Runway  layout  (similar  to  O' Hare 
International  Airport)  used  in  the  simulations 


Models  21-25  are  "shear-line"  models;  an  abrupt  wind  shift  occurs  along  a 
line.  In  models  21-24  the  winds  on  one  side  of  the  line  are  toward  the  line  at  25 
m/s  on  one  side  and  away  from  the  line  at  a  45  deg  angle  at  5  m/s  on  the  other 
side.  In  model  21  the  line  is  oriented  north-south  the  winds  toward  the  east  and 
NE.  In  model  22  the  line  is  NE-SW  with  the  winds  toward  the  SE  and  E.  For  model 
23  the  line  is  E-W  with  the  winds  toward  the  S  and  SE.  For  model  24  the  line  is 
N-S  with  the  winds  toward  the  W  and  SW.  Model  25  is  the  same  as  model  21  but  with 
an  altitude  weighting  function  that  is  zero  at  the  surface,  increases  to  unity  at 
100  m  altitude  and  decays  asymptotically  to  zero  above  100  m  altitude. 

Models  41-50  are  numerical  models  used  by  the  FAA  in  studies  of  wind- 
shear  conditions  that  were  related  to  aircraft  incidents.  They  were  furnished 
to  WPL  by  H.  Schlickenmaier  (AED-731).  These  models  are  expressed  in  a  Cartesian 
coordinate  system  centered  at  the  touchdown  point,  (x^,  y  ,  0).  The  winds  in 
the  Cartesian  coordinate  system  defined  in  section  5  are  found  by  coordinate 
transformation  and  linear  interpolation  is  used  to  obtain  values  needed  at 
points  scanned  by  the  radar. 

Models  41,  42,  43,  and  45  are  one -dimensional  models  —  vertical  profiles 
of  horizontal  wind  with  no  vertical  velocity.  Model  41  is  a  logarithmic  profile 
that  represents  the  neutral  boundary  layer.  Models  42  and  45  are  frontal  wind 
models  related  to  incidents  at  Tokyo  and  Boston.  Model  43  is  from  a  profile 
along  the  glide  path  in  a  thunderstorm. 

The  other  models  are  two-dimensional  wind  fields  which  represent  the  hori¬ 
zontal  and  vertical  winds  in  a  vertical  plane  through  the  approach  path.  Models 
44,  46,  47,  48  and  50  represent  thunderstorm  winds  and  model  49  is  for  a  frontal 
wind.  The  data  are  supplied  as  vertical  wind  profile  at  various  horizontal 
distances  from  touchdown. 

For  models  1-5  and  21-25  the  wind  speed  and  direction  are  the  same  for  all 
simulations  and  the  viewing  direction  of  the  radar  is  varied  by  selecting  different 
runways.  This  shows  how  the  radial  component  of  the  same  wind  profile  would 
appear  with  different  observation  angles.  For  models  6-12,  the  location  of  the 
downburst,  and  therefore  the  wind  speed  and  direction,  are  varied  depending  on 
tie  runway.  For  models  41-50,  the  wind  input  is  relative  to  the  touchdown 


point  so  the  wind  field  relative  to  the  fixed  coordinate  system  and  the  radar 
scanning  direction  depends  on  the  runway  selected. 


The  sign  convention  used  is  that  positive  radial  wind  is  away  from  the 
radar,  positive  glide  slope  wind  is  a  headwind  for  a  landing  aircraft,  positive 
crosswind  causes  a  landing  aircraft  to  drift  right  and  positive  vertical  wind  is 
upward. 


7.  SIMULATION  RESULTS 


A  large  number  of  simulation  runs  were  made  (27  wind  models  X3  runways) 
with  5  or  more  output  plots  for  each  case .  Samples  of  various  configurations 
are  included  in  this  report.  Other  output  plots  are  available  on  microfilm. 

Outputs  with  other  wind  models  and  runways  could  be  simulated  with  the  flexible 
computer  program  that  has  been  developed,  wind  speeds  shown  in  figures  6.1  through 
6.21  are  in  meters  per  second. 


A.  Radial  Velocity  Profiles  (Figures  6.1-6.21) 

Model  1,  a  logarithmic  one-dimensional  profile,  has  constant  wind  direction 
from  the  west  and  wind  speed  that  increases  to  30  m/s  at  500  m  altitude.  This 
is  the  actual  headwind  that  would  be  encountered  by  a  landing  aircraft  to  Runway  9, 
([9]).  Figure  6.1  shows  the  glide  slope  headwind  as  a  solid  line  and  the  radar- 
measured  radial  velocity  along  the  glide  slope  as  diamond  points.  The  crosswind 
(shown  as  a  dashed  line)  and  the  vertical  wind  (shown  as  a  dotted  line)  are  both 
zero  for  this  model.  The  radial  velocity  changes  from  toward  the  radar  to  away 
from  the  radar  as  the  azimuth  angle  passes  due  north  in  a  clockwise  scan.  Note 
that  the  radar  viewing  angle  becomes  normal  to  the  approach  path  for  runways 
9L,  9R,  4L,  and  36  (Fig.  5.2).  For  these  runways  the  sign  convention  could  be 
changed  so  that  the  radar-measured  radial  velocity  would  be  positive  for  velocities 
toward  the  radar  whenever  the  viewing  direction  was  toward  the  runway.  If  this 
were  implemented  for  Fig.  6.1,  all  the  diamond  points  would  indicate  negative 
radar-measured  radial  velocity,  with  a  zero  at  the  point  where  the  radar  viewing 
angle  was  normal  to  the  approach  path.  However,  if  there  is  crosswind,  the  radial 
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velocity  may  become  zero  at  any  azimuth  angle  or  it  may  never  be  zero  depending 
on  the  magnitude  and  direction  of  the  crosswind.  Changing  the  sign  convention 
would  then  introduce  a  discontinuity  in  the  profile,  so  this  was  not  implemented. 
The  identification  at  the  bottom  of  each  figure  shows  the  wind  model,  the  runway, 
the  coordinates  of  the  touchdown  point  and  the  coordinates  of  the  radar.  The  storm 
location  for  models  6-12  is  the  location  of  the  center  of  a  downburst.  The  last 
identifier  shows  whether  the  radar  measurements  were  deriyed  by  integrating  over 
the  pulse  volume  (INTT)  or  if  the  value  at  the  center  of  the  pulse  volume  was  used 
(INTF) .  Figures  6.2  and  6.3  show  how  the  same  wind  is  encountered  by  aircraft  and 
observed  by  radar  along  the  approaches  to  [4]  and  [14].  Crosswind  and  headwind  are 
present  on  both  paths.  A  radar  at  the  end  of  the  runway  would  measure  the  solid 
line  in  each  figure.  If  the  radar  had  also  been  used  to  measure  wind  profiles 
as  discussed  in  section  3,  the  diamond  points  could  be  predicted  from  the  profile 
measurements.  In  addition,  for  this  wind  model  and  all  other  one-dimensional 
models  the  radar  would  have  measured  the  glide  slope  wind  (the  solid  line)  for  each 
runway  if  the  radar  antenna  had  been  pointed  parallel  to  the  approach  path,  ihe 
radial  velocity  measurements  along  the  glide  slope  may  not  show  as  abrupt  a  change 
with  altitude  as  the  wind  components  (see  Fig.  6.2)  because  of  the  viewing  angle. 
Note  also  that  if  the  wind  field  had  been  constant  in  speed  and  direction  with 
altitude,  the  radar-measured  radial  velocity  would  still  show  a  change  with  height. 
Thus  the  rapid  change  from  negative  to  positive  velocity  at  low  altitude  in  Fig. 

6.1  is  primarily  a  result  of  change  in  observing  direction.  The  radial  profile,  by 
itself,  is  not  easily  interpreted.  Wind  models  3  and  5  are  also  one-dimensional 
with  constant  wind  direction.  Radial  velocity  profiles  along  the  glide  slope  are 
basically  the  same  as  those  in  Fig.  6. 1-6. 3. 

Wind  model  2  is  a  one-dimensional  profile  that  includes  wind  direction  change 
with  altitude.  Wind  is  from  the  west  at  the  surface  but  changes  to  northwest 
starting  at  30  m  altitude.  Figure  6.4  shows  the  profile  along  the  approach  to  [9]. 
The  radial  velocity  is  determined  primarily  by  the  tail  wind  below  about  30  m 
altitude  and  above  about  100  m  altitude.  From  30-100  m  the  radial  velocity  depends 
mainly  on  the  crosswind.  A  radial  profile  of  this  type  would  be  indicative  of 
shear  rather  than  just  the  change  in  viewing  direction  because  the  changes  in  slope 
of  the  radar-measured  profile  indicates  that  it  is  not  a  portion  of  a  sinusoid  as 
it  would  be  for  uniform  winds.  Below  about  100  m  altitude  the  crosswind  component 
becom<  >  significant  in  the  radar  measurement. 
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Wind  model  4,  a  one-dimensional  representation  of  cold  air  outflow  from  a 
thunderstorm,  includes  vertical  winds.  The  glide  slope  profiles  for  all  3  runways 
(Figs.  6. 5-6. 7)  are  very  different  but  each  is  clearly  indicative  of  shear  along 
the  path. 

Figure  6.8  is  representative  of  the  glide  slope  profiles  of  the  downburst 
model  (models  6  and  12).  In  simulations  with  these  models,  the  downburst  is  always 
positioned  near  the  approach  path  of  the  runway  that  is  being  scanned.  The  profiles 
are  similar  for  all  runways  because  of  the  axial  symmetry  of  the  wind  model.  The 
profile  depends  on  the  location  of  the  runway  and  downburst  relative  to  the  radar 
location . 

Figure  6.9  shows  how  the  glide  slope  profiles  change  when  the  downburst  model 
is  translated  with  a  20  m/s  wind  from  the  west.  The  radial  velocity  pattern  of  the 
downburst  (Fig.  6.8)  is  translated  by  a  portion  of  a  sinusoid.  In  Fig.  6.10  the 
translation  of  the  radial  velocity  profile  on  [9)  includes  that  portion  of  the 
sinusoid  which  is  changing  rapidly  (below  50  m  altitude)  so  the  actual  velocity 
profile  is  altered  significantly.  If  the  translation  can  be  measured  and  its 
effect  removed  from  the  radial  profile,  then  the  basic  features  of  the  downburst 
model  are  observed. 

Figure  6.11  illustrates  the  simulation  results  with  wind  models  with  an  abrupt 
wind  shift  along  a  line.  The  gradual  change  in  radial  velocity  at  upper  altitudes 
is  caused  by  the  change  in  viewing  direction.  The  velocity  changes  slowly  with 
altitude  as  compared  to  [9] ,  because  the  range  to  [4]  is  greater  so  smaller  azimuth 
angles  are  used  to  scan  the  0-500  m  height  interval.  The  abrupt  change  in  radial 
velocity  would  not  be  observed  if  the  filtering  effect  of  the  radar  resolution 
volume  had  been  incorporated  in  the  simulation. 

For  models  41-50  the  input  wind  field  is  expressed  relative  to  the  runway. 
Therefore  the  actual  wind  profiles  along  the  glide  slope  are  the  same  for  all 
runways.  The  wind  fields  are  one-  or  two-dimensional  (no  variation  in  wind  in  the 
horizontal  direction  normal  to  the  approach  path) .  Simulations  using  these  models 
show  how  the  same  wind  encountered  by  the  aircraft  appears  as  a  different  radial 
velocity  profile  for  different  viewing  directions  whereas  simulations  with  models 
1-25  show  how  the  same  wind  throughout  the  airport  appears  for  various  viewing 
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directions.  Figures  6.12-6.14  illustrate  the  simulation  results  for  model  42,  a 
one-dimensional  frontal  wind.  These  figures  and  Fig.  6. 1-6. 3  illustrate  the  dif¬ 
ferences  in  the  input  wind  models  described  above.  Note  in  Fig.  6.12  that  the 
crosswind  is  so  large  that  the  measured  radial  velocity  is  determined  primarily 
by  crosswind.  The  same  wind  model,  when  observed  on  runway  [9],  shows  crosswind  at 
low  altitudes  but  not  at  higher  altitudes. 

Figure  6.15-6.17  show  the  profiles  along  [4],  [9],  and  [14]  for  model  49,  a 
two-dimensional  frontal  wind  model.  This  is  a  numerical  model  with  relatively 
coarse  resolution  and  the  linear  interpolation  required  to  transform  the  input  data 
to  wind  components  at  any  point  makes  the  actual  profile  appear  as  line  segments 
rather  than  smooth  curves.  The  radial  velocity  profiles  are  similar  for  all  3 
runways  because  of  the  similarity  in  the  structure  of  the  headwind  and  crosswind. 

Model  50  is  a  two  dimensional  numerical  model  representing  outflow  from  a 
thunderstorm  that  was  associated  with  a  June  24,  1975  accident  at  John  F.  Kennedy 
International  Airport.  Large  vertical  velocities  are  included  but  they  will  not  be 
observed  by  Doppler  radar  at  low  elevation  angles.  Hie  velocity  profiles  are  very 
similar  to  some  of  the  downburst  profiles  obtained  with  the  analytical  model  with 
trariSlation.  The  radial  velocity  profile  shows  the  strong  gradients  of  tailwind. 

Comparison  of  the  headwind  profiles  and  the  radar-measured  radial  velocity 
profiles  shows  that,  for  the  wind  models  used  in  these  simulations,  the  radar 
can  adequately  measure  the  headwind  (within  2  m/s)  even  if  there  is  some  offset 
of  the  radar  from  the  runway.  The  radar  measurements  can  be  taken  as  headwinds 
along  the  approach  to  [14]  and  the  portion  of  the  approach  to  [9]  that  is  above 
about  100  m  altitude.  Measurements  along  the  approach  to  [4]  are  strongly  in¬ 
fluenced  by  crosswinds.  It  can  be  seen  from  Fig.  5.2  that  headwinds  along  the 
approaches  to  14R,  32L,  22R,  and  27R  could  be  measured  by  the  radar,  and  head¬ 
winds  could  be  measured  along  the  approaches  to  27L,  9L,  9R  and  4L  for  that  portion 
of  the  path  that  is  above  100-200  m. 

In  addition  to  the  profiles  of  radial  velocity  and  glide  slope  winds  as  a 
function  of  height,  the  simulation  program  also  outputs  the  profiles  as  a  function 
of  distance  from  touchdown.  Figure  6.21  is  a  sample  of  this  output;  the  profiles 
are  the  same  as  those  shown  in  Fig.  6.10. 
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Figure  6.1. 


WIND  SPEED-HEIGHT  PROFILES 

Solid  lines  show  headwinds  (+)  or 
tailwinds  (-)  encountered  by  a 
landing  aircraft  at  various  heights 
above  touchdown.  Dashed  lines  show 
crosswinds  and  dotted  lines  show 
vertical  winds  along  the  approach 
path.  Positive  crosswind  causes 
a  landing  aircraft  to  drift  to  the 
right  and  positive  vertical  wind 
is  upward.  Diamond  points  are  the 
expected  value  of  radar-measured 
radial  velocity  as  the  radar  scans 
along  the  approach  path.  Positive 
radial  velocity  is  away  from  the 
radar. 


Figure  6.2. 
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Figure  6.3. 


WIND  SPEED-HEIGHT  PROFILES 

Solid  lines  show  headwinds  (+)  or 
tailwinds  (-)  encountered  by  a 
landing  aircraft  at  various  neights 
above  touchdown.  Dashed  lines  show 
crosswinds  and  dotted  lines  show 
vertical  winds  along  the  approach 
path.  Positive  crosswind  causes 
a  landing  aircraft  to  drift  to  the 
right  and  positive  vertical  wind 
is  upward.  Diamond  points  are  the 
expected  value  of  radar-measured 
radial  velocity  as  the  radar  seems 
along  the  approach  path.  Positive 
radial  velocity  is  away  from  the 
radar. 


Figure  6.4. 


Figure  6.5. 
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WIND  SPEED-HEIGHT  PROFILES 

Solid  lines  show  headwinds  (+)  or 
tailwinds  (-)  encountered  by  a 
landing  aircraft  at  various  heights 
above  touchdown.  Dashed  lines  show 
crosswinds  and  dotted  lines  show 
vertical  winds  along  the  approach 
path.  Positive  crosswind  causes 
a  landing  aircraft  to  drift  to  the 
right  and  positive  vertical  wind 
is  upward.  Diamond  points  are  the 
expected  value  of  radar-measured 
radial  velocity  as  the  radar  scans 
along  the  approach  path.  Positive 
radial  velocity  is  away  from  the 
radar. 


Figure  6.6. 
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Figure  6.7. 


WIND  SPEED-HEIGHT  PROFILES 

Solid  lines  show  headwinds  (+)  or 
tailwinds  (-)  encountered  by  a 
landing  aircraft  at  various  heights 
above  touchdown.  Dashed  lines  show 
crosswinds  and  dotted  lines  show 
vertical  winds  along  the  approach 
path.  Positive  crosswind  causes 
a  landing  aircraft  to  drift  to  the 
right  and  positive  vertical  wind 
is  upward.  Diamond  points  cure  the 
expected  value  of  radar-measured 
radial  velocity  as  the  radar  scans 
along  the  approach  path.  Positive 
radial  velocity  is  away  from  the 
radar. 


Figure  6.8. 
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Figure  6.9. 
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WIND  SPEED-HEIGHT  PROFILES 

Solid  lines  show  headwinds  (+)  or 
tailwinds  (-)  encountered  by  a 
landing  aircraft  at  various  heights 
above  touchdown.  Dashed  lines  show 
crosswinds  and  dotted  lines  show 
vertical  winds  along  the  approach 
path.  Positive  crosswind  causes 
a  landing  aircraft  to  drift  to  the 
right  and  positive  vertical  wind 
is  upward.  Diamond  points  are  the 
expected  value  of  radar-measured 
radial  velocity  as  the  radar  scans 
along  the  approach  path.  Positive 
radial  velocity  is  away  from  the 
radar. 


Figure  6.10. 
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Figure  6.11. 


WIND  SPEED-HEIGHT  PROFILES 

Solid  lines  show  headwinds  (+)  or 
tailwinds  (-)  encountered  by  a 
landing  aircraft  at  various  heights 
above  touchdown.  Dashed  lines  show 
crosswinds  and  dotted  lines  show 
vertical  winds  along  the  approach 
path.  Positive  crosswind  causes 
a  landing  aircraft  to  drift  to  the 
right  and  positive  vertical  wind 
is  upward .  Diamond  points  are  the 
expected  value  of  radar-measured 
radial  velocity  as  the  radar  scans 
along  the  approach  path.  Positive 
radial  velocity  is  away  from  the 
radar. 


Figure  6.12. 
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Figure  6.13. 


MIND  SPEED-HEIGHT  PROFILES 

Solid  lines  show  headwinds  (+)  or 
tailwinds  (-)  encountered  by  a 
landing  aircraft  at  various  heights 
above  touchdown.  Dashed  lines  show 
crosswinds  and  dotted  lines  show 
vertical  winds  along  the  approach 
path.  Positive  crosswind  causes 
a  landing  aircraft  to  drift  to  the. 
right  and  positive  vertical  wind 
is  upward.  Diamond  points  are  the 
expected  value  of  radar-measured 
radial  velocity  as  the  radar  scans 
along  the  approach  path.  Positive 
radial  velocity  is  away  from  the 
radar . 


Figure  6.14. 
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Solid  lines  show  headwinds  (+)  or 
tailwinds  (-)  encountered  by  a 
landing  aircraft  at  various  heights 
above  touchdown.  Dashed  lines  show 
crosswinds  and  dotted  lines  show 
vertical  winds  along  the  approach 
path.  Positive  crosswind  causes 
a  landing  aircraft  to  drift  to  the 
right  and  positive  vertical  wind 
is  upward.  Diamond  points  are  the 
expected  value  of  radar-measured 
radial  velocity  as  the  radar  scans 
along  the  approach  path.  Positive 
radial  velocity  is  away  from  the 
radar. 


Figure  6.18. 


Figure  6.19. 


WIND  SPEED-HEIGHT  PROFILES 

Solid  lines  show  headwinds  (+)  or 
tailwinds  (-)  encountered  by  a 
landing  aircraft  at  various  heights 
above  touchdown.  Dashed  lines  show 
crosswinds  and  dotted  lines  show 
vertical  winds  along  the  approach 
path.  Positive  crosswind  causes 
a  landing  aircraft  to  drift  to  the 
right  and  positive  vertical  wind 
is  upward.  Diamond  points  are  the 
expected  value  of  radar-measured 
radial  velocity  as  the  radar  scans 
along  the  approach  path.  Positive 
radial  velocity  is  away  from  the 
radar . 


Figure  6.20. 
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Figure  6.21. 

Wind  speed-range  profile  for  model 
9,  Runway ;9.  Range  is  distance 
from  touchdown.  Data  are  same  as 
shown  in  Fig.  6.10. 


B.  Radial  Velocity  Fields  (Figures  6.22-6.39) 

The  three  types  of  displays  of  the  radial  velocity  fields  near  the  approach 
path  are  illustrated  in  Fig.  6.22-6.24.  Figure  6.22  shows  contours  of  constant 
radial  velocity  ("isodops")  for  a  fixed  elevation  angle  (1.5  deg)  which  inter¬ 
sects  the  approach  to  [4]  at  about  250  m  height.  The  scan  is  therefore  a  section 
of  a  VAD  scan;  however,  instead  of  the  conventional  PPI  presentation  as  in 
Fig.  4.1,  where  isodops  are  radial  lines,  the  data  presentation  is  a  B-scan 
so  the  isodops  for  a  uniform  wind  appear  as  straight  lines  normal  to  the 
azimuth  axis.  Heavy  lines  denote  contour  intervals  of  10  m/s,  light  lines 
denote  5  m/s  intervals.  This  scan  intersects  the  approach  path  at  only  one 
point.  Model  1,  shown  in  Fig.  6.22  has  large  shear  at  very  low  altitudes 
which  were  not  covered  by  the  simulation  scan  so  the  contours  are  nearly  those 
for  unil  >rm  winds. 

Figure  6.23  shows  the  isodops  on  a  surface  that  includes  the  approach  path. 
The  approach  path  is  shown  with  short  vertical  arrows;  quantization  with  azimuth 
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and  elevation  increments  of  0.1  deg  causes  the  irregularities  in  the  display.  The 
display  shows  the  isodops  on  a  surface  where  the  range  and  elevation  angle  of  all 
points  with  constant  altitude  is  the  same  as  the  range  and  elevation  angle  where 
the  scan  intersects  the  approach  path.  Since  data  cure  displayed  as  a  function  of 
azimuth,  uniform  winds  will  result  in  isodops  that  appear  as  vertical  lines.  The 
low  level  wind  shear  of  model  1  is  seen  in  Pig.  6.23  because  the  simulation  scan 
covers  the  approach  path  at  all  heights  below  500  m. 

Another  display  which  includes  the  approach  path  is  shown  in  Fig.  6.24.  In 
this  case  the  antenna  scans  along  the  approach  path  so  the  data  surface  is  a  plane 
which  includes  the  approach  path,  radar  location,  and  touchdown  point.  At  each 
azimuth  angle  where  the  radar  beam  intersects  the  approaah  path,  the  radial  velocity 
is  displayed  as  a  function  of  range.  Since  all  points  along  the  range  axis  cor¬ 
respond  to  a  fixed  antenna  position,  the  isodops  for  uniform  winds  will  be  hori¬ 
zontal  lines.  The  altitude  coordinate  is  the  altitude  on  the  approach  path,  not 
the  altitude  of  all  data  points.  The  constant  elevation  angle  scan  and  display  (as 
in  Fig.  6.22)  would  be  the  easiest  to  implement  on  an  actual  radar  but  it  does  not 
observe  the  winds  along  the  aircraft  track.  The  glide  slope  scan,  (Fig.  6.24), 
could  be  implemented  on  any  computer-controlled  radar  and  the  processing  and  dis¬ 
play  would  be  as  easy  to  implement  as  for  the  constant  elevation  scan.  The  seem 
illustrated  in  Fig.  6.23  is  not  a  practical  scan  for  an  actual  radar  since  it 
displays  only  one  range  at  any  antenna  position.  The  radial  velocity  fields  shown 
for  other  wind  models  will  therefore  be  the  glide  slope  scan  where  uniform  winds 
will  have  horizontal  isodops. 

The  glide  slope  scans  for  [4],  [9],  and  [14]  using  model  4  are  shown  in  Figs. 
6.25-6.27.  The  corresponding  radial  profiles  are  in  Figs.  6. 5-6. 7.  Although  the 
isodops  clearly  indicate  that  the  wind  is  not  uniform,  only  for  [9]  do  the  contours 
show  that  the  glide  slope  crosses  a  rapidly  changing  wind.  This  is  consistent  with 
the  radial  profiles  which  only  show  large  gradients  for  [9] . 

Simulation  results  for  downburst  models  (Figs.  6.28-6.30)  show  the  rapid 
change  in  radial  velocity  along  all  3  approach  paths.  Figure  6.28  shows  the 
contours  for  [4]  with  a  uniform  wind  from  the  west  added  to  the  downburst 
model  and  Fig.  6.29  shows  the  contours  for  [9]  with  a  uniform  wind  from  the 
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south  added  to  the  downburst.  Figure  6.30  shows  the  downburst  for  [14].  The 
isodops  are  strongly  dependent  on  the  location  of  the  downburst  relative  to 
touchdown  but  the  vertical  orientation  of  the  contours  show  that  the  wind  field 
is  not  uniform. 

The  glide  slope  scans  for  [4]  using  models  44,  49,  and  50  and  shown  in 
Fig.  6.31-6.33.  Though  these  wind  models  depict  different  meteorological 
situations  (models  44  and  50  represent  thunderstorms  and  model  49  represents 
a  frontal  passage),  the  basic  characteristics  of  the  radial  velocity  fields 
are  very  similar.  Again,  the  isodops  become  oriented  toward  the  vertical 
showing  the  changing  wind  field  along  the  approach  path.  The  same  general 
character  of  the  radial  velocity  fields  are  observed  for  the  other  runways. 
Figures  6.34-6.36  show  the  isodops  for  the  same  wind  models  on  [9]  and  Figs. 
6.37-6.39  show  those  for  [14].  These  scans  suggest  that  the  orientation  of 
the  isodops  for  a  radar  scan  along  the  glide  slope  are  a  sensitive  indicator 
of  changing  winds  along  the  approach  path. 

One  of  the  problems  with  the  display  used  in  Pigs.  6.24-6.39  is  that  the 
resulting  patterns  of  radial  velocity  are  complex  and  difficult  to  interpret, 
even  for  relatively  simple  wind  fields.  A  measurement  and  display  technique  that 
might  be  suitable  for  wind  shear  warning  is  as  follows:  First,  use  the  radar  to 
measure  vertical  profiles  of  the  horizontal  wind  with  the  VAD  scan  or  4-point 
measurement .  The  winds  are  assumed  to  vary  with  height  and  linearly  with  hori¬ 
zontal  distance.  The  radar  then  scans  along  the  approach  path  and  the  measured 
radial  velocities  are  compared  with  those  that  would  be  present  if  the  radar- 
measured  wind  profile  were  applicable.  The  difference  field  is  then  displayed. 

If  the  difference  field  is  less  them  seme  threshold  value,  the  profile  would  be 
assumed  to  apply  along  the  approach  path.  Areas  where  the  difference  field  ex¬ 
ceeded  some  threshold  could  simply  be  coded  as  a  shaded  area  on  a  display.  For 
example,  using  this  technique  only  the  2  and  3-dimensional  wind  models  would 
have  a  difference  field.  There  would  be  no  need  for  operator  interpretation  — 
the  system  would  yield  a  wind  profile  and  an  indication  of  when  that  mean  profile 
may  be  assumed  to  apply  along  a  particular  approach. 
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Figure  6.22. 

Contours  of  constant  radial  veloc¬ 
ity  ("isodops")  for  a  fixed  antenna 
elevation  angle  scan.  The  scan 
intersects  the  glide  path  at  only 
1  point.  Uniform  winds  would 
appear  as  vertical  lines. 
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Figure  6,23. 

Isodops  on  surface  which  includes 
the  approach  path,  denoted  by  short 
vertical  arrows.  Radar  range  and 
antenna  elevation  angle  at  each 
altitude  are  the  same  as  the  range 
and  elevation  angle  where  the  scan 
intersects  the  glide  slope.  Uni¬ 
form  winds  appear  as  vertical 
lines. 
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Figure  6.24. 


Contours  of  constant  radial  veloc¬ 
ity  in  a  slant  plane  which  includes 
the  glide  slope  (shown  by  short 
vertical  lines),  the  location  of 
the  radar,  and  the  touchdown  point. 
Uniform  winds  appear  as  horizontal 
lines  because  all  points  along 
the  range  axis  (at  fixed  altitude) 
are  measured  at  a  fixed  azimuth 
angle.  The  altitude  coordinate 
shows  the  altitude  on  the  approach 
path  but  not  the  altitude  of  all 
data  points.  Heavy  lines  are 
-1 

10  m  s  contour  intervals?  light 
lines  are  5ms1"  intervals. 


Figure  6.25. 


Figure  6.26. 


Contours  of  constant  radial  veloc¬ 
ity  in  a  slant  plane  which  includes 
the  glide  slope  (shown  by  short 
vertical  lines) ,  the  location  of 
the  radar,  and  the  touchdown  point. 
Uniform  winds  appear  as  horizontal 
lines  because  all  points  along 
the  range  axis  (at  fixed  altitude) 
are  measured  at  a  fixed  azimuth 
angle.  The  altitude  coordinate 
shows  the  altitude  on  the  approach 
path  but  not  the  altitude  of  all 
data  points.  Heavy  lines  are 

10  m  s  1  contour  intervals;  light 

lines  are  5ms1  intervals. 


Figure  6.27. 


Figure  6.28. 


Contours  of  constant  radial  veloc¬ 
ity  in  a  slant  plane  which  includes 
the  glide  slope  (shown  by  short 
vertical  lines) ,  the  location  of 
the  radar,  and  the  touchdown  point. 
Uniform  winds  appear  as  horizontal 
lines  because  all  points  along 
the  range  axis  (at  fixed  altitude) 
are  measured  at  a  fixed  azimuth 
angle.  The  altitude  coordinate 
shows  the  altitude  on  the  approach 
path  but  not  the  altitude  of  all 
data  points.  Heavy  lines  are 

10  m  s  1  contour  intervals;  light 

lines  are  5ms1  intervals. 


Figure  6.29. 


Figure  6.30. 


Contours  of  constant  radial  veloc- 
5  ity  in  a  slant  plane  which  includes 
the  glide  slope  (shown  by  short 
vertical  lines) ,  the  location  of 
the  radar,  and  the  touchdown  point. 
Uniform  winds  appear  as  horizontal 
lines  because  all  points  along 
the  range  axis  (at  fixed  altitude) 
are  measured  at  a  fixed  azimuth 
'  angle.  The  altitude  coordinate 

shows  the  altitude  on  the  approach 
path  but  not  the  altitude  of  all 
data  points.  Heavy  lines  are 

10  m  s  1  contour  intervals;  light 

lines  are  5  ms'*'  intervals. 


Figure  6.31. 


Figure  6.32. 


Contours  of  constant  radial  veloc¬ 
ity  in  a  slant  plane  which  includes 
the  glide  slope  (shown  by  short 
vertical  lines),  the  location  of 
the  radar,  and  the  touchdown  point. 
Uniform  winds  appear  as  horizontal 
lines  because  all  points  along 
the  range  axis  (at  fixed  altitude) 
are  measured  at  a  fixed  azimuth 
angle.  The  altitude  coordinate 
shows  the  altitude  on  the  approach 
path  but  not  the  altitude  of  all 
data  points.  Heavy  lines  are 
-1 

10  m  s  contour  intervals;  light 
lines  are  5ms1  intervals. 


Figure  6.33. 


Figure  6 . 34 


Contours  or  constant  radial  veloc¬ 
ity  in  a  slant  plane  which  includes 
the  glide  slope  (shown  by  short 
vertical  lines) ,  the  location  of 
the  radar,  and  the  touchdown  point. 
Uniform  winds  appear  as  horizontal 
lines  because  all  points  along 
the  range  axis  (at  fixed  altitude) 
are  measured  at  a  fixed  azimuth 
angle.  The  altitude  coordinate 
shows  the  altitude  on  the  approach 
path  but  not  the  altitude  of  all 
data  points.  Heavy  lines  are 
-1 

10  m  s  contour  intervals;  light 
lines  are  5ms1  intervals. 


Figure  6.35. 


Figure  6.36. 


Contours  of  constant  radial  veloc¬ 
ity  in  a  slant  plane  which  includes 
the  glide  slope  (shown  by  short 
vertical  lines),  the  location  of 
the  radar,  and  the  touchdown  point. 
Uniform  winds  appear  as  horizontal 
lines  because  all  points  along 
the  range  axis  (at  fixed  altitude) 
are  measured  at  a  fixed  azimuth 
angle.  The  altitude  coordinate 
shows  the  altitude  on  the  approach 
path  but  not  the  altitude  of  all 
data  points.  Heavy  lines  are 

10  ni  s  1  contour  intervals;  light 

lines  are  5ms1  intervals. 


Figure  6.37. 
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Figure  6.38. 


Contours  of  constant  radial  veloc¬ 
ity  in  a  slant  plane  which  includes 
the  glide  slope  (shown  by  short 
vertical  lines) ,  the  location  of 
the  radar ,  and  the  touchdown  point . 
Uniform  winds  appear  as  horizontal 
lines  because  all  points  along 
the  range  axis  (at  fixed  altitude) 
are  measured  at  a  fixed  azimuth 
angle.  The  altitude  coordinate 
shows  the  altitude  on  the  approach 
path  but  not  the  altitude  of  all 
data  points.  Heavy  lines  are 

10  m  s  1  contour  intervals;  light 

lines  are  5ms1  intervals. 


Figure  6.39. 
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8.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  major  objectives  of  the  NOAA-FAA  program  for  investigating  the  use 
of  pulse  Doppler  radar  for  wind  shear  warning  at  airporjts  were  to:  (a)  demon¬ 
strate,  with  modifications  to  existing  FAA  radars,  that  pulse  Doppler  radar 
can  be  used  to  measure  winds  in  the  optically  clear  air  as  well  as  in  preci¬ 
pitation.  (b)  Explore  the  limitations  of  Doppler  radar  as  a  wind  shear 
detector,  and  (c)  determine  how  a  single  Doppler  radar  should  be  used  for  wind 
shear  detection. 

Experimental  work  at  NAFEC  with  an  ASR-8  radar  equipped  with  a  15  foot 
diameter  parabolic  antenna  has  demonstrated  that  the  radial  wind  component 
can  be  measured  in  various  weather  conditions .  There  has  not  been  sufficient 
experimental  evidence  to  show  that  this  radar  can  measure  the  wind  in  all 
weather  conditions,  particularly  in  the  clear  air  during  the  winter  months 
when  radio  refractive  turbulence  is  decreased.  The  experimental  work  has 
been  limited  to  measuring  radial  wind  profiles;  vertical  profiles  of  the 
horizontal  wind  and  radial  velocity  fields  have  not  been  measured.  The 
existing  hardware  at  NAFEC  is  sufficient  to  perform  these  tests.  Measurements 
of  wind  profiles  can  be  made  with  the  data  system  used  for  the  wind  shear 
program,  and  measurements  of  radial  velocity  fields  (scanning  antenna)  can 
be  made  with  the  data  system  that  was  developed  for  storm  turbulence  detection. 

In  orler  to  fully  exploit  the  ASR-8  system  at  NAFEC  as  a  real-time  wind 
measurement  system  the  data  rate  needs  to  be  increased  so  that  about  32  range 
intervals  zould  be  measured  simultaneously .  This  would  require  a  Fast  Fourier 
Transform  (FFT)  array  processor  for  the  minicomputer  system.  A  color  display 
and  multi-noment  display  would  also  be  needed.  While  the  present  system  is 
sufficient  for  verifying  potential  capability  with  off-line  processing  and 
limited  reil-time  processing,  it  lacks  the  features  required  for  an  operational 
demonstration.  If  the  FAA  decides  to  continue  to  pursue  the  use  of  existing 
radars  for  weather  information,  then  real-time  capability  will  have  to  be 
added  for  ooth  wind  shear  and  storm  turbulence  work.  The  role  of  existing 
radars  depends  on  the  role  that  new  generation  weather  radars  (NEXRAD)  will 
have  for  airport  meteorology.  A  NEXRAD  radar,  if  suitably  located  and  devoted 
to  airport  surveilance,  could  provide  storm  surveilance  and  storm  turbulence 


indications  for  the  terminal  controllers.  The  radar  would  probably  be  suffi¬ 
ciently  remote  from  the  landing  area  so  the  airport  itself  would  not  be  in 
the  clutter  pattern.  This  would  preclude  its  use  for  wind  measurements  with 
useful  height  resolution,  but  it  would  be  an  important  part  of  a  wind  shear 
warning  system  because  it  would  monitor  storms  as  they  approach  the  airport. 
Another  radar  located  at  the  airport  and  solely  devoted  to  wind  measurements 
would  be  one  possible  solution  to  the  airport  weather  problem. 

Simulation  results  show  that  a  radar  with  a  150  m  range  resolution  and 
a  1.5  deg  beamwidth,  located  at  the  airport,  is  capable  of  resolving  low-level 
wind  profiles  that  are  hazardous  to  aircraft  during  landing  or  takeoff. 

Vertical  profiles  of  the  horizontal  wind  at  the  radar  site  can  be  measured 
with  an  antenna  elevation  angle  of  3  to  10  deg.  It  should  be  possible  to 
measure  the  profile  to  several  hundred  meters  altitude  in  all-weather  condi¬ 
tions,  and  much  higher  in  precipitation.  Whether  the  wind  profile  at  the 
radar  is  the  profile  that  will  be  encountered  along  the  aircraft  track  can 
be  tested  by  measuring  the  radial  velocity  profile  along  the  glide  path  and 
comparing  it  with  the  expected  profile  based  on  the  measured  vertical  profile 
of  horizontal  wind.  A  suitably  located  radar  could  measure  headwind  profiles 
along  one  or  more  approach  paths.  [If  measurement  of  headwind  along  the 
actual  approach  path  is  sufficient  for  a  wind  shear  warning  system,  an  alternative 
to  a  single  radar  devoted  to  wind  shear  would  be  to  have  several  transmitter- 
antennas  served  by  a  single  data  system  so  that  the  headwinds  could  be  monitored 
along  all  approach  and  departure  paths.  The  sensitivity  of  each  radar  would 
not  have  to  be  as  great,  the  antennas  would  be  stationary  and  such  a  system 
could  be  competitive  in  cost.  An  FM-CW  radar  could  be  considered  for  this 
application. ] 

A  single  Doppler  radar  cannot  measure  the  actual  wind  along  the  approach 
path  in  the  complex  meteorological  situations  associated  with  convective 
storms.  In  this  case  the  radar  can  only  infer  wind  shear  from  the  measurement 
of  the  radial  velocity  and  its  gradients  along  the  approach  path  and  the 
surrounding  area.  It  should  be  noted  that  vertical  winds  cannot  be  measured 
along  the  approach  path  by  a  Doppler  radar  since  the  elevation  angles  will  be  too 
low  to  sense  a  significant  component  of  vertical  motion.  Detection  of  dangerous 
vertical  winds  must  be  indirect  —  that  is,  their  presence  must  be  sensed  by  the 


perturbations  (divergence)  of  the  accompanying  horizontal  flow.  Simulations 
show  that  wind  shear  can  be  detected  by  a  scanning  Doppler  radar.  T*»e  simulations 
used  did  not  investigate  whether  every  wind  shear  could  be  detected  from  radial 
velocity  measurements  but  it  did  show  that  shear  in  specific  models  used  to 
depict  dangerous  winds  could  be  detected. 

A  shear  detection  and  display  technique  for  operational  use  is  suggested 
for  a  single  Doppler  radar  which  is  located  so  that  it  cannot  observe  headwinds 
on  all  approach  paths.  It  would  measure  wind  profiles  which  are  assumed  to  apply 
only  at  the  radar  location.  (The  temporal  changes  in  these  profiles  are  one  in¬ 
dication  of  the  lateral  extent  of  their  validity.)  The  radar  would  then  compare 
the  radial  velocity  measurements  made  by  scanning  the  area  around  the  approach 
path  with  those  expected  from  the  profile  measurements.  If  the  measurements 
agreed  the  profile  would  be  assumed  to  be  valid  along  the  approach  path.  If  they 
did  not  agree  (within  specified  limits),  then  the  wind  profile  would  not  apply 
and  the  area  could  be  shown  simply  on  a  control  display  as  a  region  of  wind  shear 
caused  by  horizontal  wind  gradients.  This  type  of  display  and  warning  system 
could  be  investigated  further  using  the  computer  programs  developed  for  the  work 
reported  here. 
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APPENDIX:  COMPUTER  SIMULATION  PROGRAM 


Hie  structural  diagram  for  the  expanded  simulation  program  is  given  in  Fig. 

12 

A.l,  and  the  associated  input-output  table  is  Table  A.l  .  The  program  has  been 
expanded  to  simulate  radar  scanning,  to  display  3-dimensional  wind  components,  and 
to  include  several  additional  wind  models.  In  addition,  seme  modifications  were 
made  in  an  attempt  to  make  the  program  somewhat  more  simple  and  more  efficient. 

For  example,  some  of  the  modules  were  combined  to  eliminate  redundancy,  some  param¬ 
eters  that  were  input  previously  (such  as  radar  beam  width  and  glide  path  angle) 
are  now  defined  as  constants,  and  some  parameters  that  were  passed  through  several 
levels  of  hierarchy  are  in  a  COMMON  block. 

In  the  fixed  beam  simulation  used  in  Phase  I,  the  origin  of  the  Cartesian 

12 

coordinate  system  was  taken  as  the  touchdown  point  .  In  the  present  study,  how¬ 
ever,  it  was  desirable  to  have  the  same  radar  scan  several  different  runways,  so 
that  a  more  natural  choice  for  origin  is  the  radar  location.  The  expanded  program 
allows  for  a  completely  arbitrary  coordinate  system,  including  the  possibility  of 
the  radar  and  runways  being  at  different  altitudes. 

The  following  parameters:  touchdown  coordinates,  XT,  YT,  ZT;  the  radar  coordi¬ 
nate,  XR,  YR,  ZR;  the  wind  reference  coordinates,  XC,  YC,  ZC;  the  runway  number, 

NRW  (for  aircraft  heading);  and  a  flag  INTP  to  select  radar  radial  velocities 
either  at  a  point  or  integrated  over  pulse  volume,  are  passed  by  the  COMMON  block, 
POSIT.  In  the  current  program  INTP  is  also  used  to  determine  if  a  volume  con¬ 
taining  the  glide  path  is  scanned  (point  velocity  option)  or  a  scan  intersecting 
the  glide  path  issued  (integrated  velocity  option) .  The  wind  reference  coordinates 
are  used  to  locate  the  center  of  the  downburst  models  or  an  arbitrary  point  on  the 
shear  line  for  those  particular  models. 

As  in  Phase  I  the  model  number  MDN  is  passed  by  COMMON/MODEL/  and  the  spe¬ 
cialized  graphics  routines  are  omitted  from  this  discussion.  Due  to  the  length  of 
the  expanded  program  it  is  impractical  to  include  a  complete  computer  program 
listing  in  this  report,  but  it  is  available  from  the  authors  upon  request. 
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STRUCTURAL  DIAGRAM 


Computer  simulation  program  structural  diagram 


Module  description 


GUSIM  -  driver  that  initializes,  controls  and  terminates  program 

execution . 

PARAMS  -  reads  in  parameters  for  each  approach  and  radar  scan. 

PROFILE  -  computes  Cartesian  coordinates  along  glide  path  for  radar 

scam  and  controls  scan  and  velocity  contour  plots. 

HCOMAX  -  determines  maximum  wind  component  for  profile  plot  scaling. 

RPLOT  -  plots  wind  component  velocity  profiles  as  a  function  of 

aircraft  range  (distance  from  touchdown) . 

ZPLOT  -  plots  wind  component  velocity  profiles  as  a  function  of 

aircraft  altitude. 

LIMITS  -  determines  limits  and  grid  spacing  for  velocity  contour 
plots . 

GPCOMP  -  computes  actual  wind  velocity  components  along  glide  path. 

RADCOR  -  transforms  Cartesian  coordinates  into  spherical  radar 

coordinates  for  radar  seem. 

RMOOMP  -  computes  radar  radial  wind  velocity  components  along  the 
glide  path. 

GPSCAN  -  computes  radial  velocity  components  for  radar  scanning 

a  volume  that  contains  the  glide  path. 

GPSCTN  -  contours  radar  radial  velocity  fields  as  a  function  of  radar 
azimuth  and  aircraft  height  and  as  a  function  of  radar 
range  and  aircraft  height;  and  plots  the  glide  path  on 
each  contour  plot. 

FESCAN  -  computes  and  contours  the  radial  velocity  field  as  a  function 

of  radar  range  and  azimuth  at  a  fixed  elevation  angle  that 
intersects  the  glide  path  at  about  5  km  from  touchdown. 

SCALE  -  scales  the  profile  plots. 

RNINT  -  computes  the  radar-measured  mean  radial  velocity  within  a 

pulse  volume  by  integrating  in  range  the  mean  radial  velocity 
in  the  beam  cross  section. 

CONTUR  -  scales  and  draws  contour  plots. 

ELINT  -  computes  the  mean  radial  velocity  in  a  beam  cross  section  by 

integrating  with  respect  to  elevation  the  mean  radial  velocity 
in  the  azimuth  beamwidth. 
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AZINT 

RADVEL 

WIND 
WIND  (N) 

TDPROF 


computes  the  mean  radial  in  the  azimuth  beamwidth. 
converts  the  spherical  coordinates  to  Cartesian  and  computes 
the  radial  component  of  u,  v,  sold  w. 
provides  decision  table  for  selecting  wind  model, 
computes  the  wind  components  u,  v,  and  w  at  the  point  (x,y,z) 
for  wind  model  N  (N=l,  ...  ,  12,  21,  ...  ,  25,  Dl,  ...  ,DX) . 
interpolates  the  numerical  models  41  through  50  to  the 
point  (x,y,z)  and  transforms  components  to  u,v,w. 


The  routines  WINDDl  through  WINDDX  are  for  models  41  through  50. 


INPUT-OUTPUT  TABLE 


INTER-  INPUT  OUTPUT 

FACE 


1 

MDN 

Model  number 

2 

GPW 

Headwind  component 

GPC 

Crosswind  component 

GPV 

Vertical  wind  component 

GRG 

Distance  from  touchdown 

GPZ 

Aircraft  altitude 

NCP 

Number  of  points  in  glide 
path  profile 

RMW 

Radar  radial  wind  velocity 
along  glide  path 

NRM 

Number  of  points  in  radar 
wind  profile 

3 

GPW, 

GPC,  GPV,  NGP,  RMW,  NRM 

WMX 

Absolute  maximum  wind 
velocity  component 

4 

WMX, 

NRM 

GPW,  GPC,  GPV,  GPR,  NGP,  RMW, 

5 

WMX, 

NRM 

GPW,  GPG,  GPV,  GPZ,  NGP,  RMW, 

6 

CSEA 

Array  of  sines  and  cosines 

EIM 

Mid-point  elevation  angle 

of  glide  path  azimuth  and 

AA 

Radar  azimuth  grid  points 

elevation  angle 

RA 

Radar  range  grid  points 

GPR 

RMN 

Minimum  radar  range 

NA 

Number  of  radar  azimuth (and 
range)  points  in  volume  scan 

7 

X,  Y, 

Z  Cartesian  coordinate  of 

GPW, 

GPC,  GPV  At  the  point 

the  point 

CSEA 

8 

X,  Y, 

Z 

RR 

Radar  range 

RMN 

REL 

Radar  elevation  angle 

DRM 

RAZ 

Radar  azimuth 

9 

RR,  REL,  RAZ,  DRM 

RMW 

At  a  point 

drr 

Range  integration  increment 

BW 

Radar  beamwidth 
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10 

RR,  REL,  RAZ,  RA,  AA,  NA 

NH  Number  of  height  increments 
in  the  radar  scan 

AG  Radar  azimuth  of  the  scan 

and  glide  path  intersection 

EG  Radar  elevation  angle  of 

the  intersection 

RG  Radar  range  of  the  inter¬ 

section 

RZ  Height  of  the  center  of 

pulse  volume  nearest  the 
glide  path 

VRA  Radial  velocity  matrix  as 
a  function  of  height  and 
azimuth 

VRR  Radial  velocity  matrix  as 
a  function  of  height  and 
range 

11 

AA,  RZ,  AG,  NA,  NM 

VR  is  VRA  or  VRR 

KT  Scan  type  indicator 

12 

RA,  ELM,  AA 

VRA  Radial  velocity  matrix  as 

NR  Number  of  range  steps  in 

a  function  of  azimuth  and 

radar  scan 

height 

13 

WMX 

YUP  Maximum  y-axis  value 

YLO  Minimum  y-axis  value 

KDY  Number  of  y-axis  grid 

divisions 

14 

WMX 

XUP  Maximum  x-axis  value 

XLO  Minimum  x-axis  value 

KDX  Number  of  x-axis  grid 

divisions 

15 

x,  y,  z 

U,V,W  Cartesian  wind  components 

16 

EE,  REL,  RAZ,  DRR,  DRM 

Mean  radial  velocity  within  the 
pulse  volume 

17 

RR,  REL,  RAZ 

Radial  velocity  at  a  point 

18 

R  is  RA  or  RR 

AZ  is  AA  or  RAZ 

REL 

Radial  velocity  at  a  point 

19 

AA,  RZ,  NA,  NH,  KT 

V  is  VRA  or  VRR  from  volume  scan 

20 

AA,  RA,  NA,  NR,  KT 

V  is  VRA  from  fixed  elevation  scan 

21 


RA,  AA,  ELM 


Radial  velocity 


Mean  radial  velocity  in  radar 
beam  cross-section 


R  varying  radar  range  within 
pulse  volume 

DRR,  RAZ,  RED,  BW _ 

R,  DRR,  RAZ  Mean  radial  velocity  with  respect 

EL  varying  radar  elevation  to  azimuth  within  pulse  volume 

angle  within  pulse  volume 

BW _ 

24  R,  AZ,  EL  Radar  range,  azimuth  Radial  velocity  at  a  point 

and  elevation  angle  varying 

within  pulse  volume _ 

X,  Y,  Z _ U,  V,  W _ 

X,  Y,  Z _ 0,  V,  W _ 

ZP  Array  of  heights  values  U,  V,  W 

UP  Matrix  of  headwind  component 
VP  Matrix  of  cross  wind  component 
WP  Matrix  of  vertical  wind  component 
RP  Array  of  horizontal  distances  along  glide  path 
NZ  Number  of  heights 
NX  Number  of  distances 
X,  Y,  Z 


